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ABSTRACT 


An analysis of conditions leading to atmospheric 
icing, design and calibration of an experimental icing 
Pacw iy, (FROST. Tunnel), and preliminary experiments on 
atmospheric icing are presented. The principal variables 
that characterize icing, the thermodynamic state of icing 
environment, and the mechanism of the accretion process are 
discussed. The experimental part is to check capabilities 
of the FROST Tunnel as an icing study facility rather than 
reveal the physics of icing or to obtain experimental 


design data. 


* 
FROST is an acronym for Fundamental Research On 


Solidification and Thawing. 
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CHAPTER iL 
INTRODUCTION 


The presence of water droplets and moisture in a 
cold environment generally leads to icing i.e. the form- 
ation of ice or snow on exposed surfaces. The Variety of 
surface shapes and sizes together with the variety of 
environmental icing conditions can be combined, basically, 
into three groups of practical significance and applicability 
associated with problems of: 

(1) Aircraft icing 

(2) Ship-board icing 

(3) Land-based structure icing 

Man has experienced icing and its dangers for 
centuries, especially in his early exploration of northern 
seas. Icing on ships, which is usually due to freezing 
of the sea spray that is blown over the vessel, may result 
in a heavy ice deposit on most of the above-water parts of 
a ship. This ice accretion increases the topside weight of 
the ship thus raising the center of gravity and reducing 
the eshi pestabilityvajdinespiteco£k sthisnevidence;eLtdwascthe 
dangen.of rairenat teyvscings that ebrought jresearcherseto 
systematic study of the atmospheric icing problems. 
Unfortunately, the results thus obtained could not have been 


applied to solve ship icing problems.. This led to the 
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situation’ that in the 1960's while aircraft were safely 
flying under all possible meteorological conditions man was 
alarmed by) the’ loss of trawlers in) northern’ waters (e.g. 
loss of the steam trawlers Lorella and Roderigo off the 
north coast of Iceland in January 1955). Naturally since 
then ship-board icing problem has received much more 
attention. 

Ice and frost coatings on land-based objects such as 
bridges, towers, buildings often threaten their structural 
integrity which again presents danger and hazard for man. 

The principal variables involved in atmospheric 
icing are: wind speed, liquid-water content, droplet size; 
and temperature of the environment. These differ signifi- 
cantly in cases of aircraft, shipboard, and land-based 
StLarccure @icing he tHighealtituderciouds have; relatively tow 
liguid=water content, small drop size, and low temperature. 
Moreover, these variables do not change drastically within 
such a cloud so that experimental simulation is easier to 
handle. Data on the range of meteorological variables 
pertinent to aircraft icing have been widely investigated. 
The data obtained indicate that a statistical relation 
exists between ligquid-water content, mean droplet diameter, 
temperature, and pressure altitude [1]. Icing conditions 
at sea, characterized by relatively high liquid-water content, 


large droplet size, and high temperatures, are very complex. 
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They exist as a result or combination of following: 

(1) Wind generated sea spray 

(2) Freezing rain or “drizzle “lor wet snow) 

(3) Supercooled sea fogs 

While cases (2) and (3) are more common in 
typical arctic waters, the case (1) often called 
“Spindrift” is the most serious cause of icing on ships 
in northern waters. The "spindrift characteristics" are 
very unstable and locally may vary rapidly in time depending 
on the wind~-speed, the heading of the vessel relative to 
the wind, the speed of vessel, etc. A detailed parametrical 
analysis [2] revealed that the number of variables pertinent 
LO eSsni pbeard) icing, 41s) much larger than the number of non— 
dimensional parameters characterizing the physical situation. 
Therefore experimentation on shipboard icing is an enormous 
task itself since any particular experiment can handle few 
parameters of limited range only. Although detailed and 
Consistent, experimental investigation has not yet been 
performed, there have been attempts to simply map dangerous 
regions where heavy icing can occur [3]. However, the data 
collected so far are unsufficient to draw any practical 
conclusions since there is no definite knowledge a priori 
under what conditions dangerous icing can exist. 

Freezing rain, drizzle, and wet snow or their 
combination are essentially responsible for heavy ice (snow) 


deposits on land-based objects and structures. Regarding 
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4 
this same data on frequency of freezing precipitations have 
been collected [4]. This data can be used for rough estimates 
of dangerous situations on various locations (e.g. for Canada 
see reference [5]) but again there have been no attempts at 
systematic study or experimental simulation of icing 


situations’ under these» conditions. 
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AN ANALYSIS OF CONDITIONS LEADING TO 


ATMOSPHERIC. ECING 


rel Preezing of a Supercooled 
Drop 


The atmosphere can be considered, depending on 
pressure, temperature, and humidity, as an air-water vapour- 
water-ice mixture. 

A water drop can exist ina "cold" environment 
(temperatures below freezing point of water) provided it is 
"Supercooled". Strictly, the expression "Supercooled drop" 
refers to a drop in the state of metastable equilibrium. 
This, however, is not common in literature and supercooling 
refers generally to the freezing point depression, regardless 
of the mechanism leading to it. Droplets of pure water, 
only a few microns in diameter can be supercooled down to 
-40°C in the laboratory [6]. At temperatures below -40°C 
droplets freeze spontaneously, but at higher temperatures 
they can freeze provided foreign particles - ice-nuclei are 
presented. Spontaneous freezing of pure water droplets or 
homogeneous nucleation is caused by small groups of water 
molecules becoming locked into ice-like configurations thus 
forming a nucleus for the ice phase. The probability of such 


molecular aggregates increases as the temperature is lowered. 
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6 
Heterogeneous nucleation, i.e. freezing of a water drop being 
infected by foreign nuclei, has been widely investiaged. 
Perhaps the most important result, which had not long been 
recognized, is that the freezing of a supercooled drop is a 
random process and a given drop has only a statistical 
probability of freezing at a particular temperature. For 
the median freezing temperature (temperature below which half 
the drops freeze) a linear relationship between the logarithm 
of the drop diameter (d) and the median freezing temperature 


(T) in Kelvin's scale has been discovered [7}. Thus: 
ln d= A- BT 


where A and B are constants for particular sample of water 
used. The effect of rate of cooling on freezing temperature 
has also been investigated. Bigg [7] found that with the rate 
OL cooling 0205°C/min the freezing temperature of a 1 mm drop 
was about 2°C higher than the corresponding value for the rate 
Orecooling0.5-C/min. 

The freezing temperature of a water drop is also 
influenced by the presence of salts in solution. It has been 
founa™ [oll, [9] that increasing concentration of various Salts 
decreases freezing temperature in proportion depending upon 


the particular salt. 
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2s can Nee Dynamrcseotaa Drop in 
Povo 


Consider a laminar flow of an air-water droplet 
Suspension past a certain body. Suppose the body is a 
circular cylinder and the droplets are of uniform size, 
uniformly distributed, and big enough so that their inertia 
forces can not be disregarded. Particles of air would follow 
curved streamlines approaching the cylinder; the water 
droplets would move in straighter paths because of the 
inertia forces, which will be balanced by imposed (drag) 
forces of resulting relative motion between the air and the 
water droplets. Thus, some of the water droplets would 
reach the cyclinder, others would by-pass it depending on 
their relative position to the centerline of the motion. 

The effectiveness of impaction, commonly called the collection 
GLb ciency, may—pe then expressed as the ratio. of the number 
of droplets that actually reached the body to the number that 
would have reached it had they continued in straight lines 
instead of following their respective trajectories. If Stokes 
law holds (small particles) the trajectories of water droplets 


can be calculated from equation of motion for a drop: 


Des — 
T™m ar = 67 da(v - = ee) 
hie 


where mis mass of the drop, Yr is the position vector of the 
drop, vy ic focal air velocity, d is the drop diameter, jl is 


dynamic viscosity of the air. The above equation in 
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dimensionless form becomes: 


mo me 
Sk ot SS any (2) 
at! at! 
2Re pd? 
where Ske = 
2 
9p_.R 


is Stokes number, Re is Reynolds number, p is density of 


water, Ps is density of air, R is cylinder radius, 


a a : , : 
eee R! 1.e€. normalized radius-vector, 
See eT 5 : 
ee ie ue being velocity of suspension far from 
fe) 
the body, 
Vo 
t' = Rr ty Ine Normalized time. 


The system of equations (2) was numerically 
integrated by Langmuir and Blodgett [10] for the case of 
potential flow past cylindersand spheres. Levin [11] obtained 
an analytic solution for potential flows near stagnation 
point. Also case when Re << 1 has been attempted analytically. 
The case of droplets of larger dimensions has been investi- 
gated by Fuks [12]. All the above solutions are of limited 
validity since authors completely disregarded both the 
existence of a boundary layer at the body's surface and the 
transportation of droplets by turbulent fluid beyond the 


separation point. 
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2.3 Accretion Process 

The previous discussion of nature of the freezing 
process of an individual water drop, its trajectory 
approaching an obstacle could, ‘at best, lead to a definition 
of initial cconditions"” namely collection ‘efficiency, 
thermodynamic state of the water droplets - upon which the 
accretion process begins. 

When the supercooled droplets impinge on the surface 
they begin to spread, are nucleated, and begin to freeze. 
During the freezing process the latent heat of fusion has 
to be removed and, essentially, the rate of freezing and 
temperature attained by the deposit is determined by a 
balance between the rate at which the heat is liberated by 
the freezing droplets and the rate at which the heat can be 
transferred to the environment by forced convection and 
evaporation. , Since the mean surface temperature of the 
deposit can not exceed 0°C there is for given conditions - 
ane? andwdrop “temperature, “airspeed, “and object size’ -"a 
Graltical Giiqui d-water*concentration for which all -the 
accreted droplets may be just frozen [13]. If the actual 
liguid-water concentration is lower than the critical value, 
all of the droplets freeze and the mean surface temperature 
Chithiardepositmisenbelows0nCcs “bitactual Arquid-water 
concentration is above the critical value, the excess of 
Watermwhrch cant not ibextrozen Gis! partially shed "off *mnto the 


stream and partially incorporated into: the ice structure. 
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10 
Thus, two principal ice growth regimes have been recognized: 
Dry growth regime (liquid-water concentration is less than the 
critical), and wet growth regime (liquid-water concentration 
uspaboven ther critical the Asral result: thet followings class iti- 


cation of ice deposits has been made [14], [15]. 


Rime (dry growth regime) is a milky granular deposit 
of ice formed by rapid freezing of supercooled water droplets 
as they impinge on objects and freeze mainly in an individual 
Manner. It is lighter - the density could be as low as 0.1 
g/em> - softer, and less transparent than glaze, discussed 
below. Factors which favour rime formation are small drop 
size, slow accretion, high degree of supercooling and rapid 
removal of latent heat of fusion. Rime can be subdivided 
into two groups; soft and hard rimes, the latter generally 


Teeerring oO e-rime deposited by denser Supercooled clouds. 


Glaze (wet growth regime) is a relatively transparent 
ice deposit formed by freezing of a film of supercooled 
water. Glaze is denser, harder, and more transparent than 
rime. Factors that favour glaze formation are large drop size, 
rapid accretion, less supercooling, and slow removal of 
latent heat of fusion. Glaze is often described as compact 
ana spongy. Compact Glaze mq’ density close to 0.9 g/em? kiss 
formed when the supercooled water droplets have time to 
spread over the surface and form continuous film before 


freezing and heat transfer between the deposit and the 
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11 
environment is such that the surface is in a "just-wet" 
CONCEELON;, ~1.e% aAt@0®¢?+and all accréted water freezes. 
spongy glaze - density 0.9 — 1.0 a/ent - is formed when 
heat transfer between the environment and the deposit is not 
sufficiently rapid to allow all of the deposited water to 
freeze;,: i.e. only’ a*fraction of it freezes immediately to 
produce a skeletal framework of ice that may retain three 
to four times its own mass of unfrozen water. 

Consider a circular cylinder of diameter R moving 
with velocity V in a supercooled cloud of uniform droplet 
Size and distribution: let w be the liquid water concentration 
of the cloud and E the collection efficiency of the 
deposition process. The rate at which heat is liberated by 
freezing droplets per unit length of the cylinder is 


Qy = 2R vw [L. + oe ies - T,) + Cc; (Ty - T,)) 


where L- is latent heat of fusion, oe and c, are specific 
heats of water and ice, respectively, Tor Tes T 4 are 
temperatures of air, freezing, and deposition, respectively. 
On the other hand, the rate at which heat is transferred to 


the environment by forced convection is 
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12 
where h is the heat transfer coefficient, k is the thermal 


conductivity of the air, 
Nu = £(Re, Pr) 


is the Nusselt number, Re and Pr being Reynolds and Prandtl 


numbers, respectively. 


en 
Oo 


The rate of heat at which heat is dissipate 


evaporation is 
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= TOD L, (o5 - 0.) Sh 


where AL is the mass transfer coefficient, L. is the latent 
heat of vaporization, po. and 6 are vapour concentration at 
the deposit surface and of the surrounding air, respectively, 


D is coefficient of diffusion, 
Si Sri Sc) 


is the Sherwood number, and Sc is the Schmidt number. Hence 


the (lumped) heat balance equation is 


2RVE w[L,+c, (T -Te) ~ c. (fT -T 


¢ gq)! = 7 (k (T3-T))NutDL, (0 -o )Sh] 


d ‘o 
provided the transfer of heat by radiation has been neglected. 
The above relation could be used e.g. to calculate the mean 


surface temperature of the deposit or the critical liquid 


water concentration (Ts = 0°C) for given conditions. The 
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13 
values thus obtained (based on rather over-simplified 
relation) could be only of informative character since no 
account has been taken on the temperature at which the 
droplets actually freeze nor the heat transfer between 
freezing droplets and the underlying surface. 

Considering the accretion on a microscopic level, 
the growth rate of the ice is a function of the, so called, 
"interface supercooling" which is the difference between the 
equilibrium temperature of the ice crystal and the temperature 
of the interface. Heat transfer between the freezing or 
frozen drop is a complicated process and depends on the 
spreading of the drop on the ice surface which is a function 
of the impact speed of the drop and the temperature of 
previously formed deposit. There have been some experiments 
made, e.g. by Hallet [16] and Macklin [17] on the growth of 
ice in bulk supercooled water, but no measurments have been 
made on the growth rates of ice within supercooled droplets 
of size occuring in natural clouds. The spreading of droplets, 
conveniently expressed by a spreading factor (ratio of radius 
after spreading to initial radius), has been investigated by 
Macklin [18] and Brownscombe [19]. Recently, Macklin and 
Payne [20] attempted to formulate a theoretical model on the 
microscopic level of a deposition process associated with dry 
regime hailstone formation. They distinguish three main 
phases during which the latent heat of fusion of a freezing 


drop is removed; the initial freezing, the subsequent freezing 
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14 
and the cooling phase. The initial freezing has been described 
as a period of time when a newly arrived droplet is nucleated 
aid elus tenperature «ices rapidly towards 0°C.. The subsequent 
heeezing.is controlled partly by heat conduction into the 
deposit and partly by forced convection and evaporation to 
the environment (radiation has been neglected). When the 
freezing process is complete the droplet is cooled by forced 
convection until the initial mean temperature of the deposit 
is regained. Although there is some overlapping of these 
three phases, the total time of the heat removal has been 
regarded as the sum of the indidivual times associated with 
each particular phase. The time of initial freezing has been 
Calculated using data for the Growth velocities in buik 
supercooled water. The time of subsequent freezing and 
cooling phases have been calculated using analytical 
formulas for non-steady state heat conduction with convection 
rates based more-or-less on overaged interface temperatures. 
AS@ee Tech the rate OF unital freezing das been round 


6 


dominant (10 ° sec), compared to the subsequent freezing 


foe sec), ana the cooling phase aie SEC.) SULpELcangly, 
however, the calculation by the above theory agreed remarkably 
with the calculation based simply on the lumped heat balance 
equation. This would rather suggest that the new model of 
Macklin and Payne can not be regarded as completely realistic, 


Since there have been a number of simplifications and 


approximations, some of them of doubtful justification. For 
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15 
example only radial thermal gradient has been considered, the 
fact that at any given instant the rate of accretion and the 
heat transfer rate by forced convection vary over the surface 
has not been taken into account, the calculation in one case 
has been based on experimental data of doubtful applicability, 


CUC% 


2.4 Conclusions 

Based on the above discussion the following can be 
drawn: 

(1) Theroetical up-to-date achievements seem to be 
completely inadequate to predict intensity, shapes and 
density of ice deposition under a variety of conditions. 

(2) Experimental data for practical (design) purposes 
obtained so far are completely insufficient to be a safe 
guidance for deSigners. 

(3) Experimental results in many cases do not seem 
to be sufficient to explain the nature of the physical 
process involved. On the other hand, in many cases where 
physical mechanism of an icing process can be predicted so 
that a physical model can be established, there are no data 
on particular processes involved to support quantitative 
analyses. 

Thus it seems that experimental research on icing 
should be emphasized. This, in fact, has been the basic 
idea leading to the FROST Tunnel design which is the subject 


of the next chapter. 
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CHAPTER. Lid 
DESIGN OF THE FROST TUNNEL 


The simulation of atmospheric icing requires a 
freezing airstream containing supercooled water droplets 
of predetermined size. To simulate all possible conditions 
that occur in natural icing clouds, the following 
parameters should be controlled within an icing study 
facitaty: 

(1) Air temperature and humidity 

(2) Air speed 

(3) Air turbulence 

(4) Air pressure 

(5) Drop size and number 

(6) Drop velocity 

(7) Drop supercooling 
Thus, an! Geimg study facility isvessentially a low-= 
turbulence wind tunnel with provisions for temperature and 
humidity control and, in addition, with a water droplets 
producing system (spraying system) located within the tunnel. 

Any wind tunnel design is limited by various factors. 
Perhaps the most important ones are funds and space available. 
In the Gasetol thetrrostorannel (Figs 228 /leand!s22) the 
strict limitation was the space; 40 x 25 x 10 ft. Having the 


over-all size of the tunnel thus predetermined it was 
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1Z 
low-level turbulence requirement (large contraction ratio) 
which resulted in relatively small cross-sectional area of 
the working section (18 x 18 inches) for the basic low-speed 
(LO-SPEED in next) design. In this regard alternative 
modifications, namely, slow- and high-speed (SLO- and HI-SPEED 
in next) have been considered (see later). 

A regular octagonal cross-sectional shape of the 
working section and the return circuit was selected. This 
is an approximation of the round cross-sectional shape which 
is considered as a least power absorbent one for a given 
cross-sectional area. Also flat sides permit fewer curved 
surfaces and thus reduce manufacturing difficulties. 

The range of wind speed in the working section 
(up to 150 mph) was chosen to cover most of the conditions 
of ship-board icing, land-based structure icing , helicopter 

Lcing’, tand fa) current problemaof arctic exploration—=1cing 
On lieeingy landing, and landedjaircratt . 

A temperature range of room temperature down to 
-40°C had been originally considered. However, since the 
first stage of the experimental program will be carried out 
at higher temperatures associated particularly with ship- 
board icing the basic LO-SPEED design was equipped with a 
refrigeration unit capable of producing temperatures down 
to =20°C. Whe unit was chosen such that |the capacity can 
be increased to attain -40°C by simply adding a second stage 


compressor to the existing first stage and roughly doubling 
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the coils capacity (surface area). 


Soe Contraction and Working 


Sections 


One of the most important design factors in obtaining 
low-level turbulence and small spatial mean velocity variations 
in the working section is a large contraction ratio of the 
contraction cone located just before the working section. 

The principal advantage of a large contraction is that the 
ratio of turbulent intensity to the mean speed will decrease 
through contraction cone as the mean speed increases. 

The effects of contraction have been investigated 

* 


Bheorericatly py Prandtl. (19335) sand Waylon (1.935)... 2b U 


is the turbulent intensity defined by 


= * * 
where u is longitudional,v andw are lateral root-mean- 


square values of velocity fluctuations, respectively, then 


* 
according to Prandtl, U would vary as 


* 
This formula, for example, predicts increase of U by the 
factor of 3.3, upproximately, Lor the contraction ratio of C= 
Tw. Tnus, in this case the ratio Of the turbuient #ntensity 


to the mean speed would be reduced by the factor 0.206. In 
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Spite of the fact that this does not seem to be a substantial 
reduction most of the existing low-turbulence wind tunnels 
Rave the Contract Omeratio initne vange of 10 4to 20. The 
reasons are: Tirstlivias Prandti’s Wiheoretical) formula has 
Limitations, and secondly) a rather indirect observation, 
that a large contraction ratio results in a low air-speed 
in the settling chamber, thus permitting the installation 
of a number of damping screens without excessive penalty for 
power absorbtion. Thus, for the basic LO-SPEED design a 
COntracciaon ratio. on sloihas been used (Fig. 3.3). 

In designing the contraction cone for a wind tunnel 
the usual design condition is that the velocity at the end 
Olietne cone, belore mt neyworkingsseccion must be fairly 
uniform. However, if the curvature of the wall along the 
flow (irection 1s too large at. certain points, loca 
velocities at these points may exceed the uniform velocity 
at the end of the’contraction cone. There are then regions 
of adverse pressure gradient and the boundary layer may 
Separate from/the wall. In the particular case of a contraction 
cone, the highest velocity is reached at the wall of the cone. 
Theretorne, 1f the velocity at the wall is made;to ancbhease 
monotonically from the beginning of the cone to the end, the 
pressure at the wall will then be decreasing monotonically 
and the danger of boundary layer separation will be avoided. 
There have been a few methods developed to solve this problem 


based, essentially, on a radial expansion of an adopted axial 
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23 
velocity distribution of afiserroes,solution of the Stokes-— 
Beltrami differential equation. The method used here is 
perhaps most general, and most suitable as a practical 
design method and was suggested by Cohen [21]. The numerical 
computations were carried out using University's IBM 360/67 
system. The results for LO- (as well as SLO-) design are 


Shown in Tab. 3.13 for HI-SPEED modirtication in’ Tab. 3.2. 


Sree iiewee tt ling (Chander 

The settling chamber consists of a frame and several 
spacers, both made of aluminium rectangular channel (see Fig. 
3.4), and thus permitting installation of as many as 8 
damping screens. The use of fine mesh screens in the settling 
chamber is to improve the flow, 1.e. reduce turbulence as 
el ese spatial (mean) velocity variation in the working 
section. As shown by Schubauer [22] the theoretical reduction 


of turbulence by n damping screens is well approximated by 


ratio; 
Az Gb ts 19 


where k is the pressure-drop coefficient of one screen. 

In our casé, i.e. using 8 Screens of 20-mesh size and 0.016 
inch wire, the above ratio would be 0.0376. This would, for 
example, mean that if we have turbulence as high as 0.5% 
Without screens (for comparison N.B.S. — U.S.A. 4 1/2 Et 
famous low-turbulence wind tunnel had turbulence of 0.265% 


without screens) we could expect turbulence reduction to 
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TABLE 3.1 


LO- & SLO-SPEED Design: 
Contraction Shape Data 
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TABLE 3,2 
HI-SPEED Design: 


Contraction Shape Data 


x (inches) 0,00 ee 6,84 10,25 13.67 


ee 


y (inches) 5.00 5,00 5.00 5.03 5.05 
x (inches) 17.09 20,51 23.93 21405 30.75 
y (inches) 5,09 5,16 5,24 5.35 5.49 
x (inches) 34.18 a, 60 44.02 Ly bly 47 85 
y (inches) 5.61 5,78 6,04 6,30 6,42 
x (inches) Bure, 5 69 58,11 61.53 64.95 
y (inches) 6.63 6,90 7.09 (a Vago 
x (inches) 68,36 1 e7e P5620 78,62 82,04 
y (inches) 7.93 8,24 8.59 8.97 one 
x (inches) 85.45 88,87 92.29 95,74 O95 3 
y (inches) 10,17 10.95 11,88 43,01 15,83 
x (inches) 102, 54 105.96 109,38 112,80 416422 
y (inches) 1707 20,02 22.06 24.09 26,08 
x (inches) 119, 64 123,08 126,50 129.92 133.34 
y (inches) 27.76 29.07 Shp isp Swab 32,10 
x (inches) 136.76 140,18 143,60 147,02 150,44 
y (inches) 33.08 a5 08 34,31 34.79 55.10 
x (inches) 153.86 157,28 160.70 164,12 167, 54 
y (inches) 35.42 Boat 35.88 35.95 36,00 
x (inches) 170.96 174,38 177,80 180,00 


eee EE LT 


y (inches) 36,00 36,00 36,00 36,00 
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value of 0.019% which is, indeed, a tremendous reduction - 
far beyond actual needs. It is common knowledge that 0.1% 
turbulence or less is sufficient for any laminar aerodynamic 
experimentation in wind tunnels and has no further effect on 
transition. However, the experimental verification of 
thes therceticalWrormula’ shows’that’ particularly for’n > 5 
the formula is less accurate and predicts a bigger reduction 


than has been observed. 


JeJ Dit Lusers 

Since the power losses in a wind tunnel vary with 
the cube of the air-speed, it is desirable to increase the 
cross-sectional area of the return passage as quickly as 
possible. The rate of the area increase is limited by the 
amount of expansion which the boundary layer will stand 
without separating from the wall. Townsend [23] has shown 
that a turbulent boundary layer can survive unlimited retardation 
without separation if the velocity varies more slowly than a 
certain power of the distance from the boundary layer origin. 
This power depends on the Reynolds number but is always more 
Posicivermcian — l/s. sibhe pociuion at which the actual power 
law index falls below - 1/3 depends on the diffuser angle and 
the thickness of the boundary layer at entry, and in the 
present state of boundary layer theory precise design rules do 
not exist. Based on experimental observations, however, a 


ditfuser of 5 to 7 degs of total included angle will yield 
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28 
satisfactory flow if relatively small area ratios are used 
(up to 3) [24]. On the other hand, there are methods 
permitting use of far larger cross-sectional area expansions, 
perhaps up to 50-60 degs of total included angle: these 
methods are based on boundary layer control schemes such as 
suction and blowing, the use of splitters, vortex generators, 
windmills, screens etc. Unfortunately, the situation regarding 
design methods or experimental data available on these methods 
is even worse. Perhaps the most promising method of boundary 
layer control by suction seems to involve too many technical 
complications. 

In designing the diffusers for the FROST Tunnel the 
problem was to enalrge the size of the working section (1.5 x 
Lt5e6t2) vEO thessizelofutheesettiing chamberalGnx Gatt.) 
within the length of the tunnel available (36 ft.). There 
were three possible locations where the diverging section 
could be used. The first one, located just after working 
Sectvoners LONettsplong rand ise@thellocatronco® the first 
diffuser using the following parameters: 

Area ratio. F208 

Total angle: 5°40' 

The second diffuser, located between the first and 
second corner is 8' - 5" long and has parameters: 

Area trations s.s 


Totaleancle.: S10! 
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The third diffuser, located just downstream of the 
fan, 15 15™ —"6ledong~and thas® parameters= 

Abeal Tati9o+es. 25 

Lotalanglesy 67s" 

The fourth, wide angle diffuser,is next (downstream) 
of the third’ ditfuser’.-° The parameters are: 

Area ratio: 1.5 

Total angler: 25 °S 0 

It is clear from the above parameters and Fig. 3.2 
that the first three diffusers are of convential design i.e. 
linearly expanding ducts with parameters within recommended 
values. In fact first two diffusers must be of conventional 
cdésigm since- heavy icing May occur there so that other design 
alternatives seem to be unacceptable. The reasons for the 
conventional design of the third diffuser have been: 

(1) Boundary layer control by suction and blowing 
would present great complications (e.g. additional cooling 
system’ for “auxiliary air")! ‘and will be costly. 

(2) It is: not known precisely how and when the use 
OE splatters, windmills, Screens ,// and vortex "generators 1s 
appropriate. 9 lt} isa matter of trial and ‘error’ for’ any 
particular situation. Mereover ;an-*spite of thefactethat 
the diffuser parameters are in "safe region," the technology 
of manufacturing process may cause limitations (e.g. "wavy 
walla whichis *quitempossible’ for ‘such ‘along: pieceof duct 


made of fibreglass). In such cases the above means could be 
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used to improve performance of the diffuser. 

The reason for the particular design of the fourth 
(wide-angle) diffuser was space limitation. However, as soon 
as there is more lab space for the FROST Tunnel it will be 
replaced by another convential diffuser of the parameters: 

RV eCalrat)O.ee, > 

TOtCALeang Leo. dO. 
which will actually be a continuation of the third diffuser 
with the same parameters. 

There are two reasons for this replacement: 

(1) A wide-angel diffuser consumes more energy 
(because of screens - discussed below) than an "equivalent" 
conventional one 

(2) For FROST Tunnel modifications (see Sec. 3.6) there 
is a requirement for additional length~space of the tunnel. 

The wide-angle diffuser devised in the first place 
is a useful experiment. As stressed before, no reliable 
design data are available presently. However, Bradshaw and 
Pankhurst of N.P.L. [25] (perhaps the most qualified experts 
on wind tunnel design), say: "A 45° equivalent cone with about 
three screens fitted in intervals has been found to maintain 
satisfactory flow with an area ratio of Lis, Coubc rus 
whether the flow in wide-angle diffusers is entirely free of 
separations, but minor regions of separated flow do not seem 
to produce unsteadiness in the working section." Thus, three 


screens of 20 wire size and 16 inches wire diameter (one at 
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the beginning, middle and end, respectively) were used. 


3.4 The Corners 

The cost and limited space available in most wind 
tunnels usually requires the design of abrupt corners with 
losses kept to the minimum by means of proper turning vanes. 
Corner vane camber lines can be designed by cascade theory 
to approximate to any given camber line in an infinite 
stream, "but in practice a circular arc camber, subtending 
an angle slightly less than 90° is adequate" - as experienced 
by Bradshaw and Pankhurst [25]. It has also been found that 
sheet metal vanes are adequate, and may even be more efficient 
than streamlined vanes since excessively thick vanes in a 
cascade cause blockage and the resulting adverse pressure 
gradients over the trailing edges of the vanes will cause 
rapid boundary-layer growth which may result in separation 
and loss of efficiency. Thus a gap/chord ratio of 1/4, a 
leading edge angle of 4-5°, and trailing edge of 0° appear 
to give the best performance [25]. In the actual design the 
above rules were followed precisely. The frame and turning 


vanes were made of aluminium sheet material welded together. 


3.5 The Fan and Refrigeration 


The calculation of power and heat losses which 
resulted in the particular choice of fan and refrigeration 


unit respectively,are presented in appendix. 


yroeds abeoeso vd banpiesb od ap ‘nant senso ‘ones x9n709 =~ 
. es ist ink ws ak aget ipdmas avin yi: od asemixoxegs ot 
paibacddue YtedmBS sts “aplosaio | ‘spigossq. a ive" morte 
boomed xsqxs kp’ * "sd qupsbs eax’ *08 pa aest Adipiie aigns ae 
g64¢ Brvet asod ale aes #1 (es) seustinns ans watebost vd 


saskoe Tis s1oft ed ative yaar le. ced reine $55 say tedam toed 
& “gt sensv aac Yievisheoks sortie eensy bentimessse meet 
suueesng steliashdeasd ss on’ brs ryadsoud SeULD obaoes 


eaus. iLiw esxey ecy! to aepbes pat ini etd 120 ajasibeap 


moidsrsace ai +i uees eA do trie. age os vovgst pysebaion biqss 

5 ,8XE to ose brodo\ gen a. eudT  . canaineame to ae brs ne 

2906 eg ty aybo ii Fea Syris' , "2 Fo sigas ‘a vaibest ae 
off, roize® Lauds ‘ott at, . tec] oonsimo Bren deed eng evig oo 7 a 
ae bas svisa « eit SY Eseivexg bowel [oa arow aoften. eves | . 
_sorencnl Bonin nlite sees mpiaals 36 some eyenmmmie | 7 


a2 

3.6 FROST Tunnel. Modifications 

As mentioned before, the SLO- and HI-SPEED modifi- 
cations of the basic LO-SPEED design have been considered 
for future experimentation. 

ihe SLO=SPEED fmodi fication (Fig: 2.5) is simply a 
rearrangement of two major (one existing) sections of the 
tunnel, namely, the contraction cone and the working section. 
Thus asoepy 6 ft. fand 11 ft. tong working, section feo be 
manufactured) will be located just after (downstream) the 
settling chamber whereas the existing contraction cone will 
be located after the working section. The principal advantage 
of this modification is that larger-scale model tests could 
be performed, however, in a lower speed range (0 - 15 mph). 

The *HI-SPEED modi ication APig. 2°25), as compared to 
the basic LO-SPEED design, differs in the contraction section, 
thevsize of the working section (l0ebyel0-inches) “and the 
speed range of up to 450 mph (estimated). This will enlarge 
the extent of experimental studies to include aircraft icing. 
For tie mocilacation, a contraction cone of ratio 52 and of 
Vength PS £2 ashore’ straignt—-walled ‘diffuser (lenghe 7 it., 
see Fig. 2.5) and a new working section would have to be 
manufactured. 
37 The Spraying System. The purpose of the spraying 
system is to produce a certain spectrum of droplet sizes of 
spatially uniform distribution. Each spray nozzle must be 


capable of sufficient water flow to satisfy maximum Piguid—we cee 
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34 
content requirement which was chosen to be 60 g/s. This 
represents liquid-water concentratiorof about 5 g/m? at 
maximum speed (150 mph) and would cover perhaps all possible 
icing situations. For illustration, liquid-water content 
“upeco ot g/m”, approximately, occurs in high altitude clouds. 
The range of drop srzes~and-theiredistribution “in.clouds and 
Sea Sprays (drop size.of l0imicrons up to b em) is so wide 
that it is impossible to cover the whole range by means of 
one nozzle design. Therefore, standard existing nozzles, 
capable of producing droplets of 10 to 200 microns, approxi- 
mately, were chosen for the first stage of experimentation. 
in? particular,ia “Spray Set-Up Nove consisting. of the 


"Fluid Nozzle 205" and "Air Nozzle 67147" designed and 


manufactured by Spraying System Co., Bellwood, Illinois, U.S.A. 


meets the above requirements. Considering the spray area 
coverage of one nozzle, an array of 5 nozzles was devised to 
cover uniformly the whole cross-sectional area of the working 
section. Moreover, enough flexibility in the design (the 
position and angle of inclination of the nozzles with respect 
to the stream can be changed) were left to overcome eventual 
non-uniformities in the spatial spray distribution. The 
design of the whole system is illustrated in Fig. 3.6. The 


air entering the air nozzle is preheated thus preventing 


freeze-out conditions at the water nozzle and the water circuit 


is equipped with a cooling unit to control the water temperature 


(drop supercooling). 
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Since the first stage of experimental studies will 
be devoted to icing resulting from presence of water droplets 
in atmosphere which exists naturally, under saturated 
conditions, no provisions were made for the humidity control. 
The humidity control, however, will be necessary if 
"hoarfrost" ice deposition (hoarfrost is a deposit of 
interlocking ice crystals formed by direct sublimation of 
water vapours on objects) is studied. Such:a control would 
consist of a fine spraying system and moisture removal system 
governed by a humidity (dew point) sensor: its design is being 


Garried OUL presently. 
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CHAPTER IV 


INSTRUMENTATION, CALIBRATION, AND 


PRELIMINARY EXPERIMENT 


A water-droplet-laden airstream is essentially 
defined when its velocity and temperature variation, the 
droplet size distribution and its spatial uniformity, and 


the liquid-water content, along with air turbulence is known. 


AsP) Velocity Measurement 


The wind speed was measured by means of standard 
pitot-static tube. A pitot-static tube reads the dynamic 
pressure, i.e. the difference between the total and static 


pressure, which by definition is 


where V is the airspeed, 0 is the density of air evaluated 
at temperature and pressure of the working section. The pressure 
difference was measured by an inclined fluid manometer. Usual 
corrections, namely, the tip and stem errors have been taken 
Lntoe account. 
For the velocity survey, the pitot-static tube was 
moved around and across the working section and the 
dynamic pressure was measured by a micromonometer. The 


variation of velocity in the working range of the tunnel has 
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been found to be less than 0.2% from the mean at 90 mph when 
6 dumping screens (of 20 mesh-size and 0.016 inches wire 


dia.) were in the settling chamber. 


4.2 Temperature Measurement 


It is difficult to measure temperature in an icing 
cloud and a satisfactory (direct) method has not yet been 
developed. Perhaps the best method available presently is 
to measure temperature in a low velocity water-and ice-free 
region and calculate temperature in the working section on 
the basis of the subsequent flow changes being expected; 
particularly, heating due to spray. This has the evident 
disadvantage that the actual test temperature could only be 
calculated when the amount of spray frozen is known, i.e. 
usually at the end of an experiment. However, these 
differences are not substantial and can be readily predicted. 
Hence the temoerature was measured in the settling chamber 
just upstream of the spray nozzles by means of 10 copper- 
constantine thermocouples; two of them being located at the 
centreline of the tunnel, the others circumferentially 
distributed about 5 inches from the tunnel walls and all at 
the same cross-sectional plane of the settling chamber. The 
thermocouples, equipped with cold junctions, were connected to 
a microvoltmeter through a switching unit. In calibration of 
the thermocouples, it was found (using a fine potentiometer 


and a fine mercury thermometer) that voltages read by the 


paiol ap ot sxpanrogmad exuendm of LEE ek od 
hood dey. ton gad ‘betidem (325216) yroioestias 8 fe buote 


eet. 


ai gisaeceng alan tsve. bosom vy anit <tc 7 
eaxiett bas-retsw Ytiooley wol 5 ak sxusebgmes exvensm. re 
‘gto noigoee on Eatxow sat nf Bie Tp TEqmSs ateluates bas soipet 
~bodosuxe prisd zapmsdo wel? insupssdve ails. 30, alend edt 
fmebive odt set cin? ..Yeuge of, eub pasteon \eiassuaisisg 
ed Yino bigoo siussieqeies tend fenton oA sods spetnavber tb 
26.4 ,awouwt ei. nosaxt? ysige 30 jraoms. od ety beatpivodeo 


secit .ravewol +tceriteges cs do Ene ped. ts-yp llevar 


Jbedoibexg yiihed: 2d aso hos isidnsigdve Jon sas, psoneze22ib 
socmetis enifgtesz ond ai bexsene® spy sine steqmed edt sorteH 
~xaqgoy OL 20 etsem yd acisson” yeage sit Ip aeettaqu deat 


enz° +s bstacot pated madt to ows :eSlomovorreds. satsinetenos — 


yilsitesitawosts exetto ste \Lennys en Yo entientnes 
#@ ils brs eliow fenabs sit i073 eptoni 2 tuods bosudbatekb 


oat .xedesdo prtittsa’ sis to emekd “Eenatsose-sagaa, Sina add 
St begooncao oxau saddkspautt Bkopy ddiw beqaiupe ,eokquocamieds. — 


a! 


39 
microvoltmeter followed very closely the standard (wire) 
calibration curve. 

The temperature variation across the stream was 
measured by a special probe consisting of five series - 
connected thermocouples. As a result, a temperature variation 
of less than 0.10°C was observed at -15°C of air temperature, 


and at 90 mph. 


4.3 Turbulence Measurement 

A constant temperature T.S.I. (Thermo System Inc.) 
hot wire anemometer (model no. 1051-2) equipped with a single hot 
filmsensor and a rms voltmeter were used to measure the 
turbulence intensity. Because of the availability of only 
a Single sensor, the turbulence intensity in the mean flow 
direction alone (u*/u) was measured. The probe was located 
in the working section at the centerline of the tunnel. 
The turbulence intensity was measured at room temperature 
(68°F) at 20, 50, 90, and 140 mph. Its value has been found 
to be within 0.03 - 0.08%, the latter value corresponding to 
140 mph. 
4.4 Droplet Size Distribution 

Measurement 

An oiled-slide technique [26] was used to determine 

droplet size distribution in clouds. The instrument (Fig. 4.1) 


consists of a glass slide "a" covered with a suitable oil which, 


while being exposed for a short time in a droplet-laden 
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41 
airstream, collects droplets and is then photographed 
through a microscope. In taking droplet samples, the most 
important factor to consider is the exposure time of the 
giled siide.in a cloud. If, it, is too. short, the. number. of 
droplets collected is insufficient to represent the variety 
of. sizes... if itis too leng, too many droplets are collected, 
thus increasing probability of coalescence which leads to 
an error in the observation. The proper exposure time, 
depending on the size, number density, and the velocity of 
impact, has been therefore determined by trial for each 
condition. The timing device of the instrument designed and 
used is based on the function of a hydraulic damper; the 
motion of the piston"B initiated by a spring 'C is controlled 
by the rate of by-passing of a liquid (oil) through a line 'D' 
connecting two sides of the piston inside the cyclinder. An 
oOiled-slide carrier connected to the piston has a circular 
opening "E" which moves axially inside a cylindrical housing 
with a similar opening. When these two openings coincide 
droplets are allowed to impinge on the oiled surface. The 
slide must be covered by a special oil which retards 
evaporation of droplets sufficiently long to have a 
photomicrograph taken. Golitzine[26] who examined a number 
of oils and lubricants found that the property of retarding 
evaporation does not depend on the viscosity of oil nor on 
the additives alone. However, he found that "Shell Spirax 


250" gave the most satisfactory results. Unfortunately, 
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Since this oil is not available on the market presently, 
alternative choices with the properties closest to Shell 
Spirax 250 were made. "Shell Spirax 140 E.P." and Shell 
Omala 85" were examined with the latter giving better results 
(as to the retarding evaporation) and thus being 
predominantly used. 

The photomicrographs were taken using Leitz-Wetzlar 
(Germany) microscope and Miranda (1.8/50 mm) camera loaded 
with Kodak Plus - X (125 ASA) film, the exposure time being 
1/8 sec. The prints were enlarged on 6 x 8 inches plates 
thus giving an over-all magnification of about 250. A 
semiautomatic particle-size analyser, developed and 
manufactured by the Research Institute of Chemical 
Equipment, Prague, Czechoslovakia, was used to determine 
size-distribution curves. The analyzer enables up to 1000 
droplets to be examined in one hour. 
4-58iThe Iiaquid=Water Content 

Measurement 

The nozzle water flow rate was used to determine 
the liquid-water content of a simulated cloud. The flow 
rate was measured by a mean calibrated flow-nozzle located 
before the header (Fig. 3.2) in the spraying system water 
line. However, since under some conditions, an unknown 
amount of water freezes on the tunnel walls and in addition 
some water may be lost by evaporation, a supplementary 


measurement based on the cloud-drop microphotography was 
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attempted. 

Ideatly, for a-givencwater flow rate and air-speed, 
i.e. for a given liquid-water content, (mean) drop size 
variations could have been obtained by changing the nozzle 
air pressure. As experienced, however, varying of the nozzle 
air pressure, the nozzle water flow rate and particularly 
the sairspeed soften led to a change of the ‘spray area coverage 
in the working section. Thus for every set of parameters, 
the spray area coverage had to be checked to get an accurate 


value of the liquid-water content. 


4.6 The Experiment 


The experiment was designed to check the ability of 
tue greost Alunnely te produce feonditi onsitin whichepeing Under 
both wet and dry regimes can be studied. 

A hollow circular cylinder; 2 inches in diameter and 
onet/S inchs walileenickness fixed yertically in sthetworking 
section, was selected as the subject of present icing study. 

For chosen conditions; wind speed of 46 mph (40 knots), 
air temperature of -10°C, and (mean) drop size of 70 microns, 
the calculation based on the lumped heat balance equation 
showed that the critical liquid-water content of the simulated 
cloud would be of the order of aici fan However, after a few 
trial runs it became apparent that the actual critical 
ligquid-water content for the above conditions would be higher; 


3g/m> approximately. Thus for run no. 1 (identical to 
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44 
"A4" [2]), intended to give spray glaze deposit (wet regime), 
a value of 4g /m?> was selected. 

After a few failures to obtain crystalline rime 

deposit (dry regime) - essentially because of non-uniformities 
of the spray spatial distribution in the working section - a 
"good" test (identicalto "A3" [2]), resulting in an intermediate 
(near critical) ice deposit, was completed when the (mean) 
drop size and the cloud liquid-water content were altered to 


30 microns and 4.5 g/m, respectively. 


4.7 Test Results and Discussion 

The test results are shown in Tab. 4.1, the front 
view photographs taken after 30 mins. of deposition are on 
Fig. 4.2 and Fig. 4.3, and for the photomicrographs of the 
Spray see Figs. 4.4 and 4.5. In analyzing the microphotographs, 
it has been found that the particle sizes are represented by 
dark circles whereas the lighter rings around them are 
results of light diffraction and dispersion. 

As expected run 1 resulted in a spongy glaze deposit 
(Fig. 4.3). The deposit is of a "mushroom-type" commonly 
experienced under over-critical conditions; the stagnation 
region of virtually 100% local collection efficiency was 
"warmed up" by water more intensively than any other region 
so that actual ice formulation went sidewise horizontally thus 
widening the frontal area of the deposit. As a result the 


mass icing efficiency [27] defined by ratio 
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TABLE 4.1 


TEST RESULTS 


Run Number 1 2 


Liquid-Water Content 
[g/om?] 4.0 ek 


Air Temperature 


ae 7 -10 
Wind Speed [mph] 46 46 
Mean Drop Diameter [microns] 185) 30 
Cofilection Ime yigertcul 90 FEY 
Efficiency Final 

[3] (estimated) 21.0 70 


Mass Icing: Efficiency [$] 109 22 
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weight of ice actually depos heediamel 2 shour 


weight of water spray that would deposit on projected surface 
S/o WOUL 


was relatively high (109%). 

Run No. 2, a compact glaze formation in and around 
the stagnation region with crystalline rime deposit at sides 
Of the cylinder (Fig. 4.4), was found to have A surprasingly 
low mass icing efficiency (12%). This fact brought up 
certain doubt about proper measurement of spray characteristics; 
in particular, the liquid-water content of the cloud. A 
Supplementary measurement revealed that the corrected value 
of the ‘cloud ligquid-water content was higher than it was 
before and even higher than that of run No. 1. In fact the 
flowmeter read the same value but the spray area coverage 
was smaller. 

Unfortunately, in both cases a check through 
microphetography was not possible since the timing device of 
the, sanpler did net function accurately. This ~worsened? “thc 
situation (the liquid water content was expected to be lower) 
since one would expect that a heavier spray would result in 
a heavier ice deposit. On the other hand, because of different 
(mean) drop size, the initial collection efficiency of run No. 
2) (702) was Jess than that of run Ne. 1 (90%). “Moreover, 
while the collection efficiency of run No. 2 was virtually 
not increasing during the experiment, the collection 


efficiency of run No. 1 was rapidly increasing and at the end 
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50 
of the run was estimated to be about 210%. It is doubtfull, 
however, whether one order-of-magnitude drop in the mass 

icing efficiency could be justified by the foregoing 
observation. Thus, it would appear that the actual liquid- 
water content of the simulated clouds was not. precisely 
known so that the present technique of its measurement was 


not adequate. 
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CHAPTER V 
CONCLUSIONS AND RECOMMENDATIONS 


pel CONCUSSIONS 

The calibration as well as preliminary tests revealed 
the following: 

(1) The FROST Tunnel is capable of working in the 
range of wind speeds and air temperatures for which it was 
designed (up to 150 mph and down to -20°C, respectively). 
Moreover, the range of wind speed could be increased (up 
to 450 mph) and larger scale models could be tested 
(modifications). 

(2) The airstream in the working section represents 
a fine uniform "jet" with spatial velocity variation less 
than 0.2%, turbulence intensity level in the mean flow 
direction in the range of 0.03 to 0.08%, and spatial 
temperature variation less than 0.10°C. 

(3) The spraying system is capable of producing a 
cloud liquid-water content of 60g/sec. or more. However, it 
is limited in range, depending on the wind speed and drop 
size. 

The cooling system built into the spray water line 
(water temperature control) did not function satisfactorily 
since the evaporative exchanger (Fig. 3.6) often froze. 


Also temperature and flow rate of the air flowing through a 
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line surrounding the water line were not high enough to 
prevent freezing-up conditions at the water nozzle. It was 
found that the water did not freeze in the water line but in 


the nozzle body itself. 


59.2 Recommendations 

Based on the discussion in ee previous section and 
Chapter IV the following alterations in the design, 
instrumentation,and methods used are suggested. 

(1) Spraying System: 

To improve the spatial uniformity of the spray 
in the working section an array of 9 nozzies instead of 5 
nozzles should be used. 

The evaporation exchangerof the water cooling 
system should be redesigneda or an alternative system (Fig. 5.1) 
could be considered. This consists of a single U-tube (one 
for each nozzle) immersed into the main air-stream so that 
the water temperature can be simply controlled by changing 
of depth of the U-tube insertion. 

To prevent "freeze-up" conditions at the water 
nozzles single electrical heaters may be installed within the 
nozzle bodies. 

(2) For particle-size and distribution measurement, 
although the oiled-slide technique gave satisfactory results, 
there is a need for another (checking) method. In recent 


years there have been a number of developments (e.g. based 
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on measurement of the transmission of infrared light through 
a cloud and other light-scattered techniques[ 28 J]: however, 
neleyrelvabiiity @eistill an doubt. There are a caw more 
established techniques such as direct cloud photography [29] 
and rotating multicylinder method [29], but these have 
several disadvantages. A determined search for other 
Suitable techniques is being carried out presently. 

(3) For liquid-water content measurement an improved 
timing device for the droplet sampler which would give a 
better adjustment of the exposure time will allow a required 
check through microphotography. However, since knowlecge of 
precise values of the liquid-water content seems to be an 
important factor an icing cloud simulations, yet another 
method would be desirable. Among the possible methods, the 
hot wire technigue - a wire heated electrically, exposed to 
the airstream, is cooled by impinging droplets; the magnitude 
of this cooling is a measure of the water content [30] - is 


used most widely and is recommended. 
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APPENDIX A 


POWER LOSS CALCULATION 
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60 
THE PRESSURE LOSS COEFFICIENT 


In the straight sections the pressure loss coefficient 


i.e. the’ pressure drop Ap in length L is: 


= Op L 
K = L K 5 
where D is hydraulic diameter, Kos “the skin <rrrctvon 


coefficient. For smooth pipes at high Re (Reynolds Number) 


Von Karman [31] gives: 
Ly ei 72 Th Re Oyen 8 


It is more suitable to work with a pressure loss coefficient 


related to the working section speed Ko defined by 


D 
Ko R a 
D 


where Do is the hydraulic diameter of the working section. 
Working Section 


Given: L 220 ft/sec 
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Gawverice i i= 6 35, 466 fi =) 325 TG, Re = 4025 ft/sec 
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CF COG, Ko = 0.00003 


In divergence sections (diffusers) both wall 
friction and expansion losses occur. The combined losses 


can be calculated [24] by 


: + S06 ean 0/2), tele a) 


or (3 Carve, 2 p.4 


where a is Givergence between opposite walls, Dy and Dd, 


are smaller and larger diameter, respectively. 


The First Diffuser 


61 


Given: To=sn0epemep. = 1s5 ft) Dp. = 2.5518, V5 = 7 95Lt/ sec 
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Serle, 2 2i be 10°, K 


Re D203 0 S45" 402 


2 
Hence: K = 0.0498, Ko = 0.0498 
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The Second Diffuser 


Given: I= 184256, DS = Zee tp D, lS te See ; V5 = 47 ft/sec 
Re Peet 02ies 10°, « = 0.01177 @ = 5° 10", 


Hence: K-= 070390, Ky = 07 0051 


The Third Diffuser 
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ne Wide Angle Ditruser 


There is not even an empirical formula for calculation 
of pressure losses one can rely upon. Assuming the devised 
diffuser "works" (if it does not, i.e. large separations 
occur, there |must be done "something" with it; e.g. using 
more dumping screens, anyhow) the pressure loss coefficient 
will be something in the range 0.3 to 0.5 or less. Lets 


assume the worst possible case, i.e. 
Ki 0.05, Ky 0.00044 
The Fan Lilet Discuser 


Given: L = 1.5 Et, Dy =o eo) cy Do SOO ty V5 = 40.5 ft/sec 
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The. Contraction Cone 


in” the "contraction cone=the losses «are “friction only. 


By’ Pope [24]? 


where Ly is the length of contraction cone. “Assuming a 


mean value for kK, 


Ky SO ik L/D, 
Given: La = 10 Vf, Dy =2owr oc, D, E—J0en Ars oes on Sea c= Doe Kot Oe OS 
Hence: Ko =2)'..02 7.310 


In the corners, both friction and rotation losses 
occur. For the devised guide vanes Pope [24] suggests this 
partly empirical relation: 


4.55 
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The Second Corner 


6 


Given: D Stee ner y. 


47 ft/sec Re = 0.95 x 10 


Hence: K 


0.1455 K 0.00662 


The Third and Fourth Corners 


i 


Given: D OSit Veeco ft/sec, Re-= 10.51 510 


Hence: K OL L55% Ko = 0.00605 


ihe Cooler 


As quocved by the manutacturer of the coils 


(Climatrol, Toronto) the required cooler has 


K_ = 0.078 
pie Rink eam 


For the screens Weighardt [32] suggests the following 


formuia for calculation of the resistence coefficient: 


-5/6 -1/3 
K = 6 Cl ew—sb) 5 Re 
where 
Peeared cs Olesu y= a an2 
ea total area = 1 


where 1 denotes the width of the square mesh,and d is wire 
diameter; Re is the Reynolds number based on the wire 


diameter. 


% 
‘* fs 
so 3% 


F : som ? i 


Dry 


is 

: i 
ee here % Bas aka edt 
Se ine 7 me oe 


di x 18.0 = BR .see\o? MEL = V.IdS = ¢ “emevtd 7 


ee ¥.U = a ,PBL.0 = 


sit to ustetostanss sat yao ——e. 2A 


.. ) 


; nr 
aeu iuslec0 brywiiipet 2h todaeein  ¥feo Susan) . 
avg.0 = a a 


—— a Sy AE. - 


lt 
“a 
pal 
| 
Lef 
0 
*} 

j 
nS 
*y 

_ 


ihyéedy lew armere ont sot 
7 7 mae, 
[SrneEpsItecn ebaBsSalearpe cht to not tadlvel ts ag afvmxot 


- . Tee 
C \ ie a\e- =f otk eh 


SA ata - a“ = i —) 
= 
ve 
eeD  ry e asian Io some Pa 
T BOTS fsiose 


‘ 


sitw at 6b bre eom exetpe ef2 Jo iets bee avid Beason -* he 
oan 7 


stiw sit ao’ Beeed zodmoa afloiyin ode ef oa rede 
7 
} = : , id 


The Settling Chamber Screens 
Guvensgd8=c07006 inehy-lk= 0705 Pinch eB ="0s062.oRpet= 414 


Hence: K 


aoe est Ke = 0.00495 


At most 10 screens can be used, hence 


K. = 0.0495 
Or Pk 


Pe wn rSt oabety socveen 
Sven wees e. ott, Gs 0708 Inch = Oo inch. eu 0 645 
Re = 1635 


Hence: K = 0.0906, Ko = 0.01180 
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Given: D 0; OA Sinch =) 0. anch,. 8 >= 0.6457 
Re = 835 


0.00505 


Hence: K = 0.1135, Ky 
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THE ENERGY RATIO 
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where Ay and Ve are cross-sectional area and speed of the 
testing section, respectively; p being the density of air. 
We have 
Given: p = 0.00234 slug/ft>, Ay = 1.86 me ve =8220 £24756 


Hence: 
AE 


6307, 4). tpi tuw,s 


ti 27a 
THE TOTAL STATIC PRESSURE LOSS 


Ap = 7 Or Vinge Sa 


Wee any eae 


2,96 aAnch. H,0 


SELECTION OF THE FAN 


The maximum volume flow through testing section is 


we25, 000) clm at, 150 "mph 
The total static pressure Grop 1s 


We SAInCch of HO 


A fan which meets the above requirements as well 
as the others - namely, safe operation under wet and icy 
conditions, and the price-is: 

Vaneaxial Direct Driven Fan - catalogue number 

vV 421 - X 42 - designed by 


Aerovent Fan Company, Inc., Pigua, Ohio, U.S.A., 
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manufactured under licence by 

Powlesland Engineering Ltd., Woodbridge, Ontario. 

The above fan, originally equipped with 20 hp motor, 
is a constant volume machine, i.e. regardless of air 
density changes it gives constant volume flow. As a result 
if the fan operates at lower temperatures the horsepower 
changes in direct proportion to the density. The lowest 
temperature the fan will be run at is - 40°C; hence the 
appropriate density ratio (density of air at -40°C to density 


au + 20°C) 1s 1.3. The required horsepower then as 


io ee 2m 


However, since there were available only 25 and 30 hp 


motors, the latter has been chosen. 
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HEAT LOSS CALCULATION 
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Heat Loss to the Environment (Q_) 
e 


a 


= ey maa OTe 
Os (TO rs) 
where 
= 2: 
A= 1050 Tt as overall outside surface area of 
the tunnel body. 
T. =-40°F is lowest operation temperature 
To = + /2°F 25 room temperature, and. H is overall 
heat transfer coefficient defined by 
= ot: 
Tan gers 
he k h 
et O 


where he, a are (convection) heat transfer coefficients 
related to convection inside or outside of the tunnel, 
respectively, t, k are the thickness and thermal conductivity 


of the insulation, respectively. 


hs in the above formula has been neglected since 
its relatively high value (forced convection) becomes 
unimportant with respect to the remaining terms. 

ho has been evaluated considering free convection 
from horizontal cylinders and an equivalent diameter of 


the tunnel, the latter defined by 
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actual value of De has been found 


dD. Seated ans ba Yaa ah oor 


Grashof No, 
Ceee = 25s et, 
Prandtl No. 


Pre sie 


Nusselt number [33] 


hence 


1/3 
Mo = “U54455 (Gr Pet 
= 40 
“i 
AY = 0.74 Btu/hr ft °F 


As the insulation material polyurethane was used w 


Ee = (0735 Btasae fer 


hence (for t = 1 inch): 


and 


B = 0.38 Brasher FCCor 


= a Btu/hr 
Qi, 43,900 u/ 
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POL ee ae ei cli: 

H = 0.256 Btu/hr £t7°F 
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Q3, = 20, S00. Brushy, 


Dissipation of Energy (4) 


This represents rather simple calculations; the 
tan as driven by -a S0 hp motor so atmost, all this energy 


couddvdissipate through viscous friction. 


loss could be 


Qq = 76,500 Btu/nr 
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Thus the maximum 


The maximum flow rate of water through spray nozzles 


60 gm/s 


is being considered which results in requirement of 
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ror the refrigeration unit capacity. 


Sensible Cooling of, the Spray 


Sensibie cooling of the water before the nozzles 
Q = 147200 Stu/hr 


SW 
"provided the water enters the tunnel" at +20°C. 
penst>lenCoo ling sou sthe bce, Deposit 
Oey Va 76100 Bruyn 


Sl 


provided the ice deposit temperature is -20°C. 


Hence the total heat for sensible cooling of the 


spray is 


Q. =/,21,400, Btu/hr 


ie: 


Finally the total heat capacity requirement for the 


refrigeration unit is 


199,000 Btu/hr 
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rr 
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As mentioned before since the first stage of 
experimental program will be carried out at higher temperatures, 
at most - 20°C, lower velocities, up to 90 mph,and also at 
lower water flow rates of the spray, up to 30 g/s,an 
alternative calculation of heat losses has been performed 
giving 
Heat loss to environment 

oF =) 36,000, Buu hr, 


Dissipation of energy through friction 


Og 23,000 Btuyhx 

Extraction, of Jatent hbeatstrom the spray 
Q1 = §34,200°Btu/nr 

Sensible cooling of the spray 
QO =  L2),900 Btu/hr 


S 


Hence the total heat loss 


Q 106,100 Btu/hr 


8185 tons of refrigeration 


Based on the above the following unit has been 


chosen: 
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(1) 5H40 Carrier Refrigeration Conpressor of 
Gating : 

Capacity 2.6.7 tons 

SST/502: -20°F 

Sei oe caer LO ee 

Diaive nos Lo 9 
where SST/502 x SCT/502 are saturated suction, or condensing 


temperature using refrigerant 502, respectively 


(2) 5F30 Carrier Water Cooled Condensor of 


Lacing: 


Total Meat reyectron: 13.2 tons 
Water flow rate: 16 gpm at -70°F 
Pressure drop: 6.9 ft. 


both designed and manufactured by 


Garyner Air voond1b1oning Co. 


Syracuse, New York, U.S.A. 


(3) | Cilimteakrol Evaporator 72) x. 72 inch. matched 
with the compressor; designed and manufactured by 
Climaccol Aas Coils Ltd. 


Oakville, Ontario, Canada. 


For the extension to -40°C (design data and ratings) 


use Carrier Application Data Bulletin, Form 5F, H = 2XA/1968. 
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